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MEAN FIELD DYNAMICS OF COLLISIONAL PROCESSES WITH 
DUPLICATION, LOSS AND COPY 
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Abstract. In this paper we introduce and discuss kinetic equations for the evolution of the 
probability distribution of the number of particles in a population subject to binary interactions. 
The microscopic binary law of interaction is assumed to be dependent on fixed-in-time random 
parameters which describe both birth and death of particles, and the migration rule. These 
assumptions lead to a Boltzmann-type equation that in the case in which the mean number of the 
population is preserved, can be fully studied, by obtaining in some case the analytic description 
of the steady profile. In all cases, however, a simpler kinetic description can be derived, by 
considering the limit of quasi-invariant interactions. This procedure allows to describe the 
evolution process in terms of a linear kinetic transport-type equation. Among the various 
processes that can be described in this way, one recognizes the Lea-Coulson model of mutation 
processes in bacteria, a variation of the original model proposed by Luria and Delbriick. 


1. Introduction 

The description of emerging collective behaviors and self-organization in multi-agent interac¬ 
tions started to gain popularity in the recent years and it represents one of the major challenges 
in contemporary mathematical modeling. In the biological context, the emergent behavior of bird 
flocks, fish schools or bacteria aggregations, among others, is a major research topic in popula¬ 
tion and behavioral biology and ecology [19, 20, 22, 23, 24, 32]. Other important examples of 
emergent behaviors describe building of tumors by cancer cells and their migration through the 
tissues [8, 9, 10, 11, 33]. Another famous example to consider in this contest is the classical Luria- 
Delbriick mutation problem [36, 38, 39]. Nonlinear statistical physics represents a powerful tool to 
describe these different biological phenomena. In particular, methods borrowed from kinetic the¬ 
ory of rarefied gases have been successfully employed to construct master equations of Boltzmann 
type, usually referred to as kinetic equations, describing the time-evolution of the number density 
of the population and, eventually, the emergence of universal behaviors through their equilibria 

|n. i ; ’<|. 

The building block of kinetic theory is represented by binary interactions, which, similarly to 
binary interactions between particles velocities in the classical kinetic theory of rarefied gases, 
describe the variation law of some selected agent characteristic, like its number. Then, the micro¬ 
scopic law of variation of the number density consequent to the (fixed-in-time) way of interaction, 
is able to capture both the time evolution and the steady profile, in presence of some conservation 
law [41, 43]. 

In this paper, we are interested in studying processes in which a huge population of N interacting 
agents can be characterized in terms of some scalar quantity assuming integer non-negative values, 
say (Vi,, Vn). The binary interactions between agents are described by the following rule. 
When two agents i and j interact, their pre-interaction values (V), Vj ) change to 

Vi Vj Vj Vi 

(i) r/ E- V "E v l = E A > + E Y w 

k —1 k =1 k -1 fc=1 

In (1), the quantities Xjk,Yik, Xjk,Yjk are independent random variables which assume non¬ 
negative integer values, the A’s with density function px and the Y’s with density function py ■ 

In [31], (Case 1), a collisional process of type (1), with Xjk € {0,1,2} and Yjk € {0,1}, has 
been introduced to describe the evolution of the gene-family abundance (the number of genes of a 
given family found in a genome) trough a minimal dynamics of duplication, loss and interspecies 
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horizontal gene transfer (HGT) According to this model, a fixed number N of species-genomes 
interact pairwise randomly, so that a given family can gain genes by interactions associated with 
HGT events. When two species interact they can exchange genes by HGT by drawing them from 
each other with Bernoulli trials of probability p k ■ In the same time, they draw from their own 
genome genes to be lost (with probability pi ) and duplicated (with probability p d ). In summary, 
in this case P{X ik = 0} = p h P{X ik = 2} = p d , P{X ik = 1} = 1 - Pl - p d , P{Y ik = 1} = 
1 — P{Yi k = 0} = ph- Note that, by taking p d + Ph = Pi , the mean number of genes in a given 
family is conserved in each collision. 

More generally, binary interactions (1) can be seen as a process in which the N agents can vary 
their given quantity Vi of some objects (the genes in Case 1) according to a collisional processes 
with duplication, loss and copy. An alternative interpretation of the process relies in considering 
populations in place of individuals, and individuals in place of objects. In this case, populations 
evolve following a classical branching process, and in addition interact with the other branching 
populations and exchange “migrants”. Note that, within this interpretation, the fc-th individual of 
the population i, interacting with population j, has Xi k children that remain in the population i 
(do not migrate) and Yj k that migrate to population j. 

An example of this second interpretation (Case 2) is a process in which N populations of 
cells interact and are subject to a mutation process. Here V :j represent the number of mutant 
cells in the j population. Each mutating cell can produce a clone in the same population with 
probability p and, independently, it can produce a clone which “migrates” to the other population 
with probability q. In other words P{Xi k = 2} = p, and P{Xi k = 1} = 1 — p and P{Y-i k = 
1} = 1 — P{Y ik = 0} = q. Here typically one consider situation in which the mean number of 
mutants grown in time. As we shall see this process of growth is similar to the linear process for 
cells mutation proposed by Lea and Coulson [36], which has its origin in a series of pioneering 
experiments proposed by Luria and Delbriick [39]. 

Our aim here is to make use of classical methods of kinetic theory to provide a kinetic description 
of the evolution in time of a multi-agent system obeying to binary interactions of type (1). Making 
use of this collisional mechanism between individuals, we introduce a bilinear Boltzmann-type 
equation which describes the behavior of the population in terms of its density f(y,t), where 
v represents the number of objects. Next, in the asymptotic procedure usually referred to as 
grazing collision limit, we obtain a simpler linear equation in divergence form, which retains many 
properties of the underlying Boltzmann equation, and in addition can be studied in detail. 

Having in mind as prototypes both the horizontal gene transfer model in [31] (Case 1), and 
the variant of Luria-Delbriick mutation model in [33, 47] (Case 2), we will split our analysis in 
two sub-cases, identified by the evolution of the mean density. Thus, we will limit our study to 
the cases in which the mean density of the population remains constant, or grows in time. In 
both cases, we can identify in a precise way the large time behavior of the solution. In particular, 
in presence of conservation of the mean number of objects (Case 1), it will be shown that the 
solution density converges to a steady state profile, which depends heavily from the microscopic 
interactions. 

In more details, in Section 2 we will briefly introduce the kinetic description of our collisional 
models, coupling them with some direct physical consequences. This allows the interested reader to 
take an exhaustive view of the modeling assumptions, by comparing both the nonlinear Boltzmann- 
type equation with its linear asymptotics, named quasi-invariant limit. Also, the main differences 
between Cases 1 and 2 are here underlined. Next, we collect in Section 3 the main results on 
existence and uniqueness of solutions to the mathematical models, by computing additionally 
some of the relevant mean quantities. Here, convergence to equilibrium for collisions of type (1), 
which imply conservation of the mean, value is studied in full details. Section 4 will deal with the 
large-time behaviour of the solution in the case in which collisions of type (1) imply the growth 
of the mean value. In the case of growth, by suitably scaling the solution with respect to its 
time-dependent mean value we will show that the scaled solution converges towards a fixed steady 
profile as time tends towards infinity. Last, Section 5 will describe and justify from a mathematical 
point of view the quasi-invariant asymptotics procedure leading from the Boltzmann-type equation 
to the linear one. Some technical results are for the sake of readability postponed to the Appendix. 
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2. The kinetic equation and its quasi-invariant collision limit 


2.1. The kinetic equation. Performing the usual mean-field approximation when the number 
of particles (colliding entities) goes to infinity, the evolution in time of the number density can 
be quantitatively described by a bilinear Boltzmann-type equation. In this equation, the time 
variation of the density ft(v) = f(v, t), with v G N and t > 0, follows from a balance between gain 
and loss terms, that, for the given number v, take into account all the interactions of type (1) 
which end up with the number v (gain term) as well as all the interactions which, starting from 
the number v, lose this value after interaction (loss term). This Boltzmann-type equation reads 
[43] 

( 2 ) ft{v) = Q + {ft,ft){v)-f t (v ). 

In (2), Q + is the gain collision operator defined, for any pair of densities / and g on N, by 


(3) 


w u 2 

Q + (/, 9) (v) = Prob ( E + E = v ) ’ 


*=l 


i=l 


where V\. V 2 . X 1 , X 2 , ■. ., Yj, Y 2 ... are stochastically independent, Vi has density /, V 2 has density 
g, the Xj’s have the same law of a random variable X with density px and the Y)’s have the same 
law of a random variable Y with density py. 

Equation (2) is coupled with an initial condition fo(v), which is here assumed to be a probability 
density on N. Considering that the variable v can only assume values in N, equation (2) can be 
fruitfully rewritten in terms of probability generating functions (p.g.f.) 


f t{z ):=^z v f t (v) 2 €[ 0 , 1 ], 

•u>0 


A standard computation shows that ft{z) satisfies the (simpler) nonlinear equation 

Tu/t(X) = ft{Px(z))ft(p Y {z)) - ft(z) z G [0,1], t > 0, 

(4) v- 

fo(z) =^2z v fo{ v ) 

v>0 

where 

Px{z) = E[^ A '] = Y, z m Px{m) and p Y (z) = E[z r ] = £ z m p Y {m). 

rri> 0 m> 1 

Equation (4) is the analogous of the Fourier transformed Boltzmann equation for maxwell pseudo¬ 
molecules, introduced in kinetic theory of rarefied gases by Bobylev [13]. 

To avoid trivial situations, in what follows we shall assume that both P{X = 0} ^ 1 and 
P{Y = 0} ^ 1. 

Existence and uniqueness of solution of equation (4) can be proven in a standard way (cf. 
Section 3.1 for details). 

In this paper, we will only consider random variables X and Y such that 

(5) E[X r + Y r ] < +00 

for some r > 1. In presence of condition (5) for some integer r > 1, and assuming that the initial 
density fo(v) has bounded moments up to order r, one can easily reckon explicit expressions for 
the evolution of the moments of ft(v). In particular, the mean J Z v v ft( v ) of f t evolves 

according to 

(6) M 1 (f t ) = M 1 (f 0 )e a ' t , 

where or := E[(X + Y)] — 1. Indeed, recall that for any p.g.f. p (of a probability p on N) 

Mi(p) ■=Y vp ^ = dzp ( z ) l*=i> 
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where d z p{z)\ z —\ := lim z ^\- d z p{z) < +oo if and only if Y^ V V P( V ) < +oo. By taking the 
derivative with respect to z on both sides of (4), and evaluating the resulting equation in z = 1, 
one obtains 

±M 1 (f t ) = a 1 M 1 (f t ), 

which gives (6). Analogous computations can be done to evaluate higher-order cumulant functions 
of the density ft,. An explicit expression of the variance of ft, is presented in Proposition 6. 

Note that the condition 

(7) E[X + Y] = 1 

implies that the mean remains constant in time, i.e. Mi(/ t ) = m 0 for every t > 0, as in Case 1. 

2.2. The quasi-invariant collision limit. The bilinear Boltzmann-like equation (2), fruitfully 
written in the form (4), is the starting point to obtain, in a well-established asymptotic procedure, 
simpler models that are reminiscent of the binary collision rules, and, while maintaining most 
of the properties of the nonlinear kinetic model, result to be linear. This kind of asymptotic 
procedure is close to the so-called grazing collision limit for the Boltzmann equation [48], and has 
been widely used in kinetic theory to obtain Fokker-Planck type equations to describe, among 
others, the cooling in granular gases [28, 42], wealth distribution in a multi-agent society [21], 
price formation [16, 17], and opinion formation [26, 46]. 

This asymptotics procedure is based on the following assumptions. Given a small positive 
parameter e, assume that the random variables X and Y are defined by 

(8) X= m X + (1 - m) and Y= m Y, 

where X , Y , rj i, rj 2 are independent random variables, and, for some constants b±, b 2 >0 such that 
bit < 1, i = 1,2, it holds 

P{ m = 1} = 1 - P{r) % = 0} = he, i = 1,2. 

Note that by construction the laws of X and Y depend on e, so that X = X e and Y = Y e . 
Nevertheless, for the shake of notational simplicity, the dependence in e has been omitted. Since 
e is assumed to be small, the post-collision quantities remain equal to the pre-collision 

ones (Vi,V 2 ) with high probability. Indeed, with high probability, X e is equal to one and Y e to 
zero. In order to evaluate the mean value of ft.(v) in correspondence to the collision defined by 

(8) , observe that 

E[X e + Y e ] - 1 = e[6i(E[A] - 1) +6 2 E[T[] =: ech. 

Hence, (6) takes the form 

(9) Mi(ft) = d z f t ( 1) = m 0 e aiet . 

This shows that, in presence of the small parameter e, to observe the same variation of the mean 
value corresponding to e = 1, it is necessary to wait a longer time r, given by r = t/e. 

By virtue of (8) one obtains 

px(z) = z + eb^p^z) - z) and p Y (z) = 1 + e6 2 (Pv( 2 ) ~ !)• 

Therefore, if ft(z) = ft.,e(z) is a solution of (4) corresponding to random variables X and Y defined 
in (8), and the initial density fo has finite mean mo, by expanding f t (px{z)) and ftfpyiz)) in 
Taylor’s series around z (around 1, respectively) gives 

ft(px(z)) = ft(z) + ebi(j)x(z) - z)d z f t (z) + ebiR he (t,z), 

and 

Mpy(z)) = ft( 1 ) + eb 2 (pY{z) - l)d z f t {l) + eb 2 R 2 ,e{t, z), 

where t and R 2 e denote the remainders. Since by definition f t ( 1) = 1, recalling also (9), the 
collision term in (4) takes the form 

ft(Px{z))f t (p Y (z)) - f t (z) =e b 2 f t (z){pY(z)-l)M 1 (f t ) + b 1 (px(z)-z)d z ft{z) +eR 3>e (t,z), 
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where i? 3e (f, z) denotes the new remainder term. In view of the previous remark on the evolution 
in time of the mean value, let us set r = t/e and g T ,e{z) = . f T / e {z ). Then, g T ,e{z) satisfies the 
equation 

(10) d T g T ^{z) = g T ,e{z)b 2 {Py( z ) - 1 ) M 1 (g T ,e) + h(px(z) - z)d z g T>e (z) + i? 3ie (r/e, z). 

Let us assume that f? 3>e -) 0 as e -> 0. Since Mi(< 7 T;e ) = m 0 e aiT , letting e —>■ 0, and defining 
g T (z) = lim e ^ 0 g T , e (z) shows that g T satisfies the equation 

drg T {z) = m 0 e aiT b 2 (pY(z) - 1 )g T (z) + h(px(z) - z)d z g T (z), 

(11) 

9o{z) = fo(z). 

While the actual derivation of equation (11) is largely formal, every step can be made rigorous. 
We postpone a precise statement of the limit procedure and its detailed proof to Section 5 (cf. 
Proposition 10). 

Here it is worth noticing that from the previous equation it follows that the mean of g T (') is 
moe“ lT . Indeed recall that y(r) := d z g T (z )\ z =i = M-\ {g T {-)), so that deriving with respect to z 
equation (11) and evaluating the resulting equation for z = 1 one gets 

y(j ) = 6i(E[X] - l)y(r) + b 2 m 0 e^ T E[Y]. 

This implies 

y (r) = m 0 e air + e blE[rlr (y(0) - m 0 ), 
and the claim follows since the initial condition is y( 0) = mo- 

2.3. The Lea-Coulson equation for mutation processes. Consider the special case in which 
P{X = 2} = 1 — P{X = 1} = p and P{Y = 1} = 1 — P{Y = 0} = q. With this choice 
Px(z) — z = pz(z — 1) and Py(z) — 1 = q(z — 1). By setting /?i = pbi + qb 2 , f3 2 = b\p and 
p = b 2 qmo, equation (11) takes the form 

(12) d t g(t, z) = (z- 1) {P 2 zd z g(t , z) + pe 0lt g(t, z)} . 

Equation (12) is referred to in the literature as the Luria-Delbriick model. The model deals with 
the estimation of mutation rates, and has its origin in a series of classic experiments pioneered 
by Luria and Delbriick [39], which were at the basis of the construction of a mathematical model 
able to estimate them. The original model proposed by Luria and Delbriick assumed deterministic 
growth of mutant cells, which seemed a too stringent assumption to allow for efficient extraction 
of reliable information about mutation rates from experimental data. This shortcoming of the 
model of Luria and Delbriick was some year later remedied by a slightly different mathematical 
formulation proposed by Lea and Coulson [36], who adopted the Yule stochastic birth process to 
mimic the growth of mutant cells. In the ensuing decades the Lea-Coulson formulation occupied 
a so prominent place in the study of mutation rates that the Lea-Coulson formulation is now 
commonly referred to as the Luria-Delbriick model [50]. 

The classical linear procedure of growth of both normal and mutant cells, that give rise to 
equation (12) has been merged into the framework of kinetic theory in [33]. There, the mutation 
problem was described in terms of linear interactions, and the solution to the underlying linear 
Boltzmann equation was shown to converge towards the solution to the Lea-Coulson model. A 
precise statement of this convergence result has been subsequently done in [47]. 

While equation (12) coincides with the mathematical formulation of Lea and Coulson [50], it is 
important to remark that the meaning of the various (eventually time dependent) coefficients in 
(12) assume a different meaning depending on the original model from which this equation comes 
from. 

In the original formulation, equation (12) is recovered through the following assumptions. The 
process starts at time t = 0 with one normal cell and no mutants. Normal cells are assumed to 
grow deterministically at a constant rate /3i- Therefore the number of normal cells at time t > 0 
is N(t) = Mutants grow deterministically at a constant rate f3 2 . If a mutant is generated 
by a normal cell at time s > 0, then the clone spawned by this mutant will be of size for 

any t > s. Mutations occur randomly at a rate proportional to N(t). If p denotes the per-cell 
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per-unit-time mutation rate, then the standard assumption is that mutations occur in accordance 
with a Poisson process having the intensity function 

(13) v(t) = nN(t) = 

Consequently, the expected number of mutations occurring in the time interval [0, t) is 

u(s)ds = j- (e 0lt - 1) . 

The formulation of Lea-Coulson differs from the previous Luria-Delbriick original formulation in 
that mutant cell growth is described by a stochastic birth process. If a mutant is generated by 
a normal cell at time s > 0, then at any given time t > s the size of the clone spawned by that 
mutant will have the same distribution as Y(t — s), where Y (s) is a Yule process having birth rate 
and satisfying Y(0) = 1. In this way if W(t) denotes the number of mutants existing at time 
t, then 

M(t) 

(14) W(t) = W - n), if M(t) > 1, 

i—1 

while 

W(t) = 0 if M(t) = 0 . 

Here, t, are the random times at which mutations occur, M(t) stands for the mutation process 
which is a Poisson process with intensity function i/(t) given in (13), and (T)(s)i>i is a sequence 
of independent Yule processes. In particular, the probability generating function of W ( t ) satisfies 

(12) [50]. 

Unlike the previous approach, the Lea-Coulson equation (12) follows here by linearizing a 
nonlinear equation of Boltzmann type based on classical binary interactions, and consequently 
from a nonlinear mechanism, in which the number of cells is subject only to duplication (the 
variable X) and mutation (the migration variable Y). In particular, the initial condition is here 
different, and the coefficients /?i and f3 2 are no more independent. In addition, the exponential 
in time term in front of the probability generating function g is in our derivation generated by 
the evolution of the mean value. Nevertheless, from a mathematical point of view, this derivation 
could share a new light on the large-time behaviour of the solution to (12), which is now connected 
with the solution of a nonlinear Boltzmann-type equation, that, unlike the linear model, possesses 
solutions of self-similar type. This new connection will be studied in more details in a companion 
paper. 

2.4. Steady states. When E[Y + Y] = 1, while E[X r +Y r ] < +oo for some r > 1, one can prove 
that, for any given constant too, the kinetic model (2) possesses a steady state, namely a density 
satisfying 

(15) /oo =0 + (/oo,/oo), 

of mean mo- Moreover, any solution f t with initial density /o of mean value mo converges (point- 
wise) to the corresponding stationary solution whenever Yh v vr foi v ) < + 00 (cf. Theorem 3). 
In Section 3.1 we shall show that ft converges to exponentially fast in time with respect to 
suitable Fourier metrics (cf. Theorem 2). 

If E[Y 2 + Y 2 ] < +oo, one can recover an explicit formula for the second moment of 

/oo. Indeed by combining Lemma 19 in Appendix with identity (15) one gets 

M 2 (/oo) = M 2 (Q + {foo,foo)) = M 2 (f 00 )[E[X] 2 +E[Y] 2 ] 

+ 2E[A]E[Y]M 1 (/ 00 ) 2 + [Var(X) + Var^M^). 

Since M 1 (/ oc ) = m 0 , it follows that 

(Var(X) + Var(Y))m 0 
n/ooJ 1-E[A] 2 -E[Y] 2 ' 

While in principle one can recursively determine every integer moment of the steady state 
an explicit general expression for the p.g.f /A, is not available, with the exception of very simple 
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situations. An explicit example is obtained by fixing P{X = 1} = 1 — P{X = 0} = p and 
P{Y = 1} = 1 — P{Y = 0} = q with p + q = 1. In this case the steady state coincides with a 
Poisson distribution. To see this, recall that the p.g.f. of a Poisson distribution of mean mo is 
given by foo(z) = exp(rrio(z — 1))- In the special case described above, px(z) = zp + 1 — p and 
Py(z) = z( 1 — p) +p. Thus 

Q + (/oo, foo)(z) = exp(m 0 (p Y (z) - 1)) exp (m 0 (p Y (z) - 1)) = exp(mo(z - 1)) = foo(z)- 

One of the advantages of the quasi-invariant limit approximation is related to the possibility to 
obtain an explicit expression for the steady states. Indeed, if &i(E[X] — 1) + 6 2 E[F] = 0, the 
stationary equation reads 


(16) m 0 b 2 {f>Y{z) - l)<?oo(z) + h{Px( z ) ~ Z )d z g 00 {z) =0 z G [0,1] 


and gives the explicit solution 


(17) 



(1 ~P Y (s)) 
( Px ( s ) - s ) 



In Proposition 11 we will present an interesting probabilistic interpretation of (16). Other aspects 
related to the steady states will be included there. 


2.5. Steady states for horizontal gene transfer processes. Horizontal gene transfer (HGT) is 
believed to be the dominant component of genome innovation for bacteria. Recent estimates show 
that at least 32% of the genes in prokaryotes have been horizontally transferred [34, 35]. Today 
the study of the HGT can be tackled using a growing amount of genomic data. Nevertheless the 
study of mathematical models for HGT is still at the beginning. Among the few mathematical 
studies on the subject we mention the infinitely many genes model [7]. This model uses the 
coalescent theory in order to describe the underlying phylogenetic tree. A different approach has 
been considered in [31], which focuses on the description of a collisional mechanism which describes 
the time evolution of the gene-family abundance (the number of genes of a given family found in 
a genome). As already mentioned in the introduction, in the model a fixed (very large) number of 
species genomes interact pairwise randomly. When two species interact they can exchange genes 
by HGT by drawing them from each other with Bernoulli trials of probability Ph- In the same 
time, they draw from their own genome genes to be lost (with probability p{) and duplicated 
(with probability pd)- In the model different families are considered subject to a stochastically 
independent evolution. Although this simplifying assumption is probably unrealistic, the effective 
rates ( pd,Ph.,Pi ) can vary from family to family, giving rise to the observed diversity between 
families. An interesting result, which will be briefly described below, is that the model predicts 
an increased dispersion in family abundance as HGT becomes less relevant. 

In [31] it is assumed that each gene can be transferred in a single copy, that corresponds to 
assume that Y can take only the values {0,1}. In what follows we present a description of the 
steady states (in the quasi-invariant approximation) for a slightly more general situation, i.e. both 
X and Y take values 0,1, 2. For the sake of simplicity assume that b\ = b 2 = 1 and write 

P{X = k}=p k , P{Y = k}=q k k = 0,1,2. 

Since we need E[A + Y] = 1, we impose that 

Pi+qi+ 2(p 2 + < 72 ) = 1- 


With these choices it is immediate to conclude that 

(1 -Py(s)) = gi + g 2 + sq 2 
{Px( s )-s) Po-sp 2 

Recalling that (for /? ^ 0) 


a + bs b a/3 — ab 

=-p - z] + 


log ({V^)' 
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when p 0 ,p 2 > 0, (17) becomes 

„ . , mo ^ (z _i) / 1 -P 2 /P 0 \ ((91+92)P2+P092) 

9 oo(z)=e °»2' H-- - - = 

VI -P 2 /P 0 ZJ 

Note that E[X] < 1 implies both that p 0 > 0 and p 2 < po, so that 1 — P 2 /P 0 > 0. This proves that 

9oo{y) = P n 0 (z)p N {z) 

where p n 0 (z) = e m ° p ^^ z ^ is the p.g.f. of a Poisson random variable IIo of mean mo 9 2 /p 2 
and pn(z) is the p.g.f. of a Negative Binomial random variable of parameter p = P 2 /P 0 an d 
r = ^((<71 + q 2 )P 2 + Po< 72 )- Recall that a Negative binomial random variable of parameter (p, r) 

p 2 

has density 

p N (k) = ^ ^ (1 - p) r p k k = 0, 1 ,... 

In other words, in this case, the steady state is the density of the random variable 

Vyo = IIo + IV, 

where n 0 and N are independent. Since a negative binomial distribution can be represented as a 
compound Poisson distribution, i.e. 

lb 

N = J2 Li ’ 

i —1 

where IR is a Poisson distribution of mean ?Tog(l/(l — p)) and the L, ! s are independent and 
identically distributed random variables, each one with logarithmic distribution 

PL,(k) = fc[ 0 g(l)(l _ p)) ('•'=1.2....), 

one gets the representation 

ni 

loo = iio +Li. 

i—l 

If now we set g 2 = 0, we obtain 

9oo{z) =Pn{z), 

where TV is a Negative Binomial random variable of parameter p = P 2 /P 0 an d r = m-o 9 i/P 2 , he. 
of mean mo and variance TOoPo !Q\ ■ Note that this is the model introduced in [31] and briefly 
discussed above. In this case po = pi is the probability of loss of a gene, p 2 = Pd > 0 is the 
probability of duplication and q\ = ph is the probability of HGT. It is interesting to remark that 
in this case the resulting steady distribution depends only on the initial mean mo and on the ratio 
Ph/Pd- Moreover, the dispersion measured by the index V ariance / mean = 1 +pa/ph, increases 
as the probability ph of HGT decreases. The limit case is obtained when Pd = P 2 = 0. 

Indeed, if p 2 = g 2 = 0, direct computations show that ffoo(y) = exp {mo (z — 1)}, namely that 
the steady states are given by Poisson distributions. In other words Poisson distributions are 
recovered in the case in which pd = 0 , that is when no duplications occur. 

Finally, although not relevant for the HGT model, let us consider the case in which p 2 = 0, 
while q 2 > 0. Here 

Poo (y) = exp / — ^z 2 + {qi +q 2 )z - (qi + |p 2 ) }, 

l Po L 2 2 J J 

or, alternatively 

fate) = exp fe ( g )-t]}, 

where 

Pr(z) = q* 2 z 2 + q* z , q* := q 2 /{2q x + 3 q 2 ), q{ = 1 - q* 2 . 
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Since Pr(z) is the p.g.f. of a random variable R taking values 1, 2 with probability q\, q £, it follows 
that in this case the steady state is the density of the compound Poisson random variable 

n 2 

Voo 

i—l 

where II 2 , Ri, R 2 , ■ ■ ■ are independent, II 2 is distributed according to a Poisson distribution of 
mean 2 too/[po(2<?i + 3^)] and the i?;S have the same law of R. 


3. The kinetic equation: solutions, moments, steady states 

The remaining of the paper will be devoted to detailed proofs of various results about existence, 
uniqueness and asymptotic behavior of both the kinetic Boltzmann-type equation (4), and its 
quasi-invariant limit (11). 

3.1. Existence and uniqueness of solutions. The unique solution of (4) can be written in a 
semi-explicit form by resorting to the so-called Wild series [49, 43]. 

Proposition 1. Let fo{v) be a probability density. Then, the initial value problem for equation 
(4), with initial condition /o, has a unique global solution. For any given t > 0, the unique solution 
ft can be written as 

ft(z) = e -t ^(1 - e~ 1 ) n q n (z), 

n> 0 

where qo{z) = fo(z), and, for any n > 1, 

^ n —1 

Qn(z) := -J2^(P x ( Z ))^n-l-j{pY{z)). 
j—0 

Moreover, ft remains a p.g.f. for any t > 0. 

Proof. The fact that f t (z) is a solution can be checked directly. To prove uniqueness, let ft and g t 
be two solutions departing from the same initial value. Note that in place of (2) one can consider 
the same equation in integral form 

(18) f t (v) - fo(v) = [ [Q + (/s, fs){v) — f a (v)]ds v G N. 

Jo 

Define u(t) := sup z£ r 01 i \ft(z) — 9 t{z)\ < 2. Using (18) one gets u(t) < u(0) + 2 f* u(s)ds. Hence, 
Gronwall inequality implies u(t) < u(0)e 2t . Since u(0) = 0 it follows that u(t) = 0. Finally, by 
induction one shows immediately that q n is a p.g.f. for every n > 0. Hence ft is a p.g.f.. □ 

An immediate consequence of the previous result is that the unique solution to equation (2) 
can be expressed as 

(19) Mv) = e- t J2d-e- t ) n qn(v), 

n> 0 

where qo(v) = fo(v) and, for any n > 1, 


1 

n 


n— 1 

'y ', Q + {dji Qn-l-j)(v). 

j=0 


( 20 ) 


Qn{v) ■■ 
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3.2. Steady states and convergence to equilibrium. For any r > 0, let us consider the 
following metric between probability distributions over non-negative integers 


d r {f,9) 


sup 

ze( o,i) 


\m-g(z)\ 

ll-*l r 


This is a natural adaptation of the well-known Fourier distance introduced in [29]. It is clear that 
convergence with respect to the metric d r yields convergence of p.g.f.’s and, hence, point-wise 
convergence of densities. 

In dealing with general probability measures / and g on [0, +oo), it is also useful to introduce 
the weighted Laplace transform metric 


where 


4 (/,<?) := sup 1/(0 g(01 

«>o |?| 


/(0 = / e * v f(dv) and <?(£) = 


- / 


g{dv). 


It is easy to see that d r and d* are (topologically) equivalent (cf. Lemma 18). 
We remark that, in terms of the Laplace transform 


M£) = Y e * V M V )> 

V 


the Boltzmann-type equation (2) takes the form 
( 21 ) = £> 0 . 
Now, /o(C) = Yl v e ~^ v fo{ v ) is the Laplace transform of the initial density, and 
kx(0 = log(E[e _?A [) and k Y (0 = log(E[e _?y ]) 


are the cumulant functions of X and Y. 

Given a density / on N (or more generally a probability measure / on R + ), M r (f) will denote 
the r-moment of /, expressed by y)„, v r f[v) (or more generally by f v r f(dv)). To denote the 
variance of a random variable Z we will write Var(Z). Moreover, if / is a density (or a probability 
measure) Var(f) will denote the variance of a random variable Z with law /, that is Var(f) = 
M 2 (/) — Finally, for r > 1, let us set 

(22) A r := E[X] r + E[F] r and a r = E[X] r + E[F] r - 1. 


Theorem 2 (Contraction in d r ). Assume that E[X r + Y r ] < +oo for some r £ (1,2]. Let 
/ t (1) and f ( t 2) be two solutions of (2) with initial conditions /q 1 ' and f^ such that Mi(f^) = 
Mi (/q 2 ' ) ) = too- Then for every t > 0 


(23) ^(/ t (1) ,/ t (2) )<rf r (/ 0 (1) ,/ 0 (2) )e Q '- t and d* r C/f\ /f>) /f)e“- t . 

Proof. Assume that rf r (/o i/o ) < +oo, otherwise there is nothing to be proved. Recall that 
Px(z) = E[z- Y ] = E m >o z m p x {m) and p Y {z) = E[z Y ] = J2 m >o z m p Y {m). It is easy to check (cf. 
Thm. XI.1.1 in [27]) that 

(24) 1 - p x {z) = (1 -z)^2 z m P x (rn), 

m> 0 


where Px(m) = P(X > m) = 1 — 'fZ'-Lo Pj ■ The same result holds for p Y . This yields that, for 
every z £ (0,1), one can write 

l ~ p x {z) = Y z m Px(m) < Y Px(m) = E[X], 
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and |1 — p Y {z )|/(1 — z) < E[Y], Hence, for every r > 1, 

sup ll 7i Px{ l W + sup |1 ~ M ! )r =E[Xr+E[y] r = A r . 


26(0,1) 


1~Z[ 


26(0,1) 


1-zY- 


Now recall that, by Proposition 1, 


n> 0 


where q^\z) = fo \ z ) an( ^’ f° r an y n > U 


n—1 


^O) = = “ E APaWRE-jGMu))- 


i=o 


Using the bound \q) % \z)\ < 1 for every z £ (0,1) and i = 1, 2, we obtain 


< - E SU P 
n ^ 26 ( 0 , 1 ) 


kf^MuO^EiGM*)) - q?\px{z))q { 7 ) _ 1 _Ap Y {z))\ 


\l-z\ 


(25) 


< 1 y- 1 gup kf' \i>x{z)) - qf\px{z)) |1 - px(z )I 
~~ n 2e(0 P i) 


1 51 sup 

n ^*6(0,1) 


|l-PA'(2)| r |l--| r 

^-VjtM*)) - ^^(M*))! |1 ~Py{z)\ T 


|1 -py(z)\ r 


1-Z’ 


A 


n— 1 




j=0 


Hence, if d r (f^\ /q 2 '*) = d r {qtp ,q7P) < +oo then d r {q7\qi?) < +oo for every n and by Lemma 


14 




Hence, using also (62), 

d,(/ f (1) ,/ t (1) ) < E e_ ‘( 1 - e-'rdriqgKqW) 


n> 0 


< <wP,/P) E *-'o - e -T r( ^(t + i) = ‘Wo'U'V'- 1 . 


n> 0 


To recover the second inequality in (23), recall that by (67), for every r > 1, 

A* := sup kx ^ 1 '' + sup = E[Xf + E PT- 

£>o Kl £>o |?| 

Using this identity, we can repeat the previous part of proof to conclude. 


□ 


Theorem 3 (Convergence to steady states). Let E[X + Y] = 1, and E[X r + Y r ] < +oo for 
some r £ (1,2]. Then for every mo > 0, there exists a unique density foo, with mean too and finite 
moment of order r satisfying foo = Q + (foo, foo)- Moreover, provided M r (fo) < +oo, any solution 
ft departing from the initial density fo with mean mo, converges in d r -metric to the corresponding 
stationary solution foo, and 

dr(ftJoo)<d r (fo,foo)e- Mt . 

Proof. Let us observe that A r < 1 for every r > 1 whenever E[X + Y] = 1. Let fo = S mo and 
set fn+i = Q + (fn , fn ) for every n > 0. Since E[X + Y] = 1, one gets M\{f n ) = Too for every n. 
Moreover, by Lemma 19 in the Appendix and the fact that Mi(/ n ) = mo we have 

M r (f n + 1) < A + A r M r (f n ) 
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where A = mo(E[A r ] + E[F r ]) + mQ(2E[A']E[F]) r / 2 . Iterating on n yields 

n , 

(26) M r (f n+1 ) <A^K + < A” +1 < + m °' 

i =0 r 

Arguing as in (25) one gets 

dr(fn+k> fn) A d r (/n+fc— 1 5 fn— 1)A r 5 


and hence 

d r (fn+k,fn ) < dr(/fc,/o)A?. 


By Lemma 18 and (26) 

sup d r (fk, /o) < supmax{2 r+1 , J B r c r (M r (/ fc ) + M r (/ 0 ))} < C. 

fc fc 


Hence, the sequence {/ n }n>o satisfies the Cauchy condition d r (f n , f m ) < e for every n,m > 
N(e). Since Mi(/ n ) is bounded, by Prohorov theorem one obtains that /„ is weakly compact 
(tight). Combining tightness with the Cauchy condition one concludes that there is a probability 
distribution /(*, such that f n —► /oo (pointwise). Since 

fn+l(z) = f n (Px(z))f n (p Y {z)), 

taking n -> +oo, one obtains f^z) = f 0 o(px{z))f 00 (p Y (z)). The property M r .(/ (X) ) < +oo follows 
from the fact that f n —► /oo and sup„ M r (f n ) < +oo (cf. Thm. 25.11 in [12]). Uniqueness of /oo 
(in the class of distribution with finite r-moment) follows by Lemma 18 and Theorem 2. □ 


3.3. Evolution of moments. We conclude this Section by giving explicit expressions for the 
evolution of the mean and the variance of the solution f t . 

Proposition 4 (Evolution of the mean). Let f t be the unique solution of equation (4) with 
initial condition fo- If mo := Mi(/o) < +oo, then at any subsequent time t > 0 

(27) Mi{f t ) = m 0 e ait . 


Proof First of all note that 

(28) ^vQ + (fo, fo){v) = E [X + Y]m 0 = Aim 0 . 


Hence, using (20) one gets 

By Lemma 14 
(29) 

Thus, by (19) and (62), 


Ax 


Ml (q n ) = 1 V M^qj). 

71 L ' 


Mi(q n ) = mo 


i =o 

T(Ai + n ) 
r(Ai)T(n + 1) ’ 


Mi(ft) = e~* X)(l - e~ t ) n Mi(q n ) = moeT* £(1 - e"*)" r( ^,+ ^ 


71>0 


n> 0 


□ 


Proposition 5. Let E[A r + Y r ] < +oo for some r G (1, 2]. If f t is the unique solution of equation 
(4) with initial condition fo such that M\(fo) = mo and M r (fo) < +oo, then t H > M r (f t ) is a 
continuous function and 

M r (f t ) < M r (f 0 )e olrt + e art H(t). 

Here H(t) is given by 

H(t ) = [ [A r m 0 e (ai -“ r)s + (2E[X}E[Y}) r/2 m r 0 e^- a ^ s ]ds. 

Jo 
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Proof. Let us start by proving that 

(30) t ^ M r {f t ) = e -t 5^(1 - e~ t ) n M r (q n ) 

n> 0 

is a continuous function on the set {t > 0 : M r (f t ) < +oo}, which turns out to be an open interval. 
Identity in equation (30) is immediate by (19) and Fubini theorem. Clearly, the series may diverge 
for some t < +oo. We claim that there is T = T(fo, X,Y) > 0 such that M r (f t ) < +oo for every 
t £ [0, T\. Using (20) and Lemma 19 one can check that 

1 n— 1 . n— 1 j-y n— 1 

M r (q n ) = - 53 M r{Q + (Qj,Qn-i-j)) = 5Z M Mi) + — 53[ Ml fe) + 

3=0 j -0 j—0 

for a suitable B = B(X,Y,r) < +oo. Now (29) and well known asymptotics for the Gamma 
function give 

n —1 

M r (q n ) < -53 M r (qj) + Bm 13 , 
n j=o 

for suitable Bi and /3. Since Lemma 15 yields that M r (q n ) < K n+1 , (19) gives 
M r (f t ) = e-‘ 53(1 - e-^Mriqu) < e"‘ 53(1 - e"*)"^ 1 . 

n>0 n>0 

This shows that / = {t > 0 : M r (f t ) < +oo} is non-empty. Let us show that I is open. Note 
that t i-A e t M r (f t ) is non-decreasing. Set so := sup{t : M r (f t ) < +oo}. If so < +oo and 
so £ I, one can consider the function t <—> f So +t as a solution of (2) with initial condition f So 
such that M r {f So ) < +oo. By the previous argument, it must be that M r (f t+e ) < +oo for any 
e £ [0, T(f So ,X , Y)], that gives a contradiction. The continuity on I follows by observing that, for 
z = ( 1 - e ~ t ) with t £ I, £n> o z n M r (q n ) is a convergent power series. 

Finally let us prove that, if so := sup{t : M r (f t ) < +oo} < +oo then lim t _ >s - M r (f t ) = +oo. 
If So < + 00 , the previous argument show that So ^ I, that is (1 — Zq) z o M r (q n ) = +oo 

for zo = (1 — e _s °). Then, Abel theorem yields that lim._ >z -(l — z)^2 n>0 z n M r {q n ) = (1 — 
zo) En>o z oM r (q n ) = +oo. 

Now, by multiplying equation (2) by e 4 and integrating on [0, t) we obtain 

e*/t(«) = fo(v) + f e s Q + (f s ,f s ){v)ds. 

Jo 

Hence, Fubini Theorem gives 

e t M r {f t ) = M r {fo)+ [ e s M r (Q+(f s ,f s ))ds. 

Jo 

Now note that, by Lemma 19 in the Appendix 

M r (Q + (f s , f s )) < A r (Mi(/ 8 ) + M,.(/,)) + H r Mi(/ s ) r , 
where B r = (2E[X]E[Y]) r / 2 and A r = E[A] r + E[Y] r . Since M\[f a ) = moe aiS , it follows 
M r (Q+(f s , f s )) < A r m 0 e“ 18 + B r m r 0 e ra ' s + A r M r (f s ). 

Denote m r (t) = e 4 A/ r (/t). We get 

m r (t) < A r f m r (s)ds + m r ( 0) + f e s [A r moe aiS + B r moe raiS ]ds. 

Jo Jo 

Finally, thanks to Gronwall inequality which can be applied by virtue of the first part of the proof 
(see Lemma 21), it holds 

m r (t) < m r (0)e Art + f e Ar ^~ s ^e s [A r moe aiS + B r m r 0 e raiS ]ds, 

Jo 


and sq = +oo. 


□ 



14 


FEDERICO BASSETTI, GIUSEPPE TOSCANI 


Proposition 6 (Evolution of the variance). Let Var(fo) < +oo and assume that 
(31) E[X 2 + F 2 ] < + 00 . 

Then 


(32) 


where 


Var(ft) = M 2 (f 0 )e a2t - m%e 2ait + /3m 0 


L Cy,i — 0^2 


-I{ai ^ a 2 } + te a2t I{ai = a 2 } 


+ 27771Q 


g 2ai t _ e a 2 t 


L 2 aq — OL 2 


-I{ 2 ai ^ a 2 } + te a2t l{2a\ = a 2 } 


f3 :=Var(X)+ Var(Y) and 7 := E[-X]E[y]. 


Proof. Recall that, by Proposition 5, 1 1 —>• M 2 (ft) is a continuous function. Using formula (69) in 
Lemma 19 we obtain 

(33) M 2 (Q+(f s ,f s )) = M 2 (f s )[ E[Af + E[Y] 2 ] + M 1 (f s )f3 + 2 7 M 1 (/ i ) 2 . 

Combining (18), (27) and (33) one gets 

(34) M 2 (/ t ) - M 2 (/ 0 ) = [ ( a 2 M 2 (f s ) + c(s))ds, 

Jo 

with 

c(s) = mo ^/3e“ lS + 2jmoe 2aiS ^ . 

Solving the equation one obtains 

M 2 (ft) = e“ 2t (m 2 (/ 0 ) + J* c{s)e- a2S ds ), 

which implies 

Var(f t ) = e a2t (M 2 (f 0 ) + J* c(s)e- a2S ds) - mge 2ait . 

Simple computations then give (32). □ 

Note that from E[A + y] = 1 and E[A 2 + Y 2 ] < + 00 , it follows cx 2 < 0 = aq and 2y = — a 2 . 
Then (32) implies 

(35) r«r(/,) = (V ° r (- Y) + + [Vor(A) - A2°] e -M<, 

for every t > 0. 


+ 2 I 


«2 


4. The case of increasing mean 

The asymptotic behavior of the solution to the Boltzmann equation (2) when the mean value 
is preserved in time has been fully described in Theorem 3. In this case, in fact, existence of a 
steady solution together with its main properties follows. 

A completely different situation arises when aq = E[A +Y] — 1 > 0, as in Case 2. Now the mean 
value diverges to + 00 , and the usual way to extract information about the large-time behavior of 
the solution [42] is to scale the solution of equation (4), or equivalently (21), with respect to the 
mean, in such a way to maintain the mean constant. This allows to look for non-trivial asymptotics 
profiles of 

•? 




(*)' 

where 

m{t) ■.= Mi(ft) = m 0 e ait . 

It is easy to verify that the Laplace transform h of the scaled density satisfies the new equation 


(36) ^M«) + Mf)+ «.«!'.,«) = /.,(-»(( ) t *(Jn)) i >(-"'(t)C( s |g))- 
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Note that if /W and /^ are two solutions of (2), with initial conditions and /g 2 \ while 
h W and h^ are the respective scaled solutions, then 

drM X) A 2) ) < <(/t ( 1 ) ,/fV(i)“ r = <(/ t (1) jfVo r e- trQl > 

and (23) yields that 

In view of the previous bound, it appears natural to assume that 


(37) (ct r /r — or) < 0 for some r > 1. 

Condition (37) is reminiscent of the analogous condition for the existence of self-similar solutions of 
Kac like kinetic models, see [2, 14, 42]. Indeed, as we shall see, there is a precise connection between 
the asymptotic profile of the scaled solution of (2) and the asymptotic profile of the solution of 
a particular Kac like equation. It is important to remark that, since a r = E[A] r + E[Y] r — 1, 
condition (37) is equivalent to ask that r H > a r /r is decreasing in r = 1, i.e. ( a r /r )' < 0 for r = 1. 
Computing (a r /r)' one obtains that (37) is equivalent to 

(38) E[X](log(E[X]) - 1) + E[Y](log(E[Y]) - 1) - 1 < 0. 

Clearly the previous relation is true if E[Y] < 1 and E[A] < 1, while is clearly false if one of the 
two expectations is bigger than e. In Figure 1 a numerical evaluation of the region where (37) is 
verified and E[X] + E[Y] > 1 is reported. Additional information on the function r H > a r /r — or 
can be found in Lemma 3.10 in [42]. 



Figure 1. The interior of the curve (colored) represents the region in the plane 
(E[A],E[Y]) where (37) is verified and E[A] +E[Y] > 1. 


Recalling that here m(t) —»• +oo, using (66), one can write 

-m(t)k x {€/m(t)) = £E[X] + 0(£ r m(f) 1-r ), 


and hence, 

lim ( — m(t)kx{Ji/rn(t) \ =^E[A']. 

t —>- + 00 \ / 

The same result holds for Y . In this way one can see that (at least formally) h t converges to the 
solution /loo(C) = /[o,+oo) e-^’hooidv) of 

MO + ai^i^oo(0 = /ioo(PW)^oo(eE[Y]). 


( 39 ) 
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For a precise statement see next Proposition 9. Now following [2, 6, 14], (39) is equivalent to the 
integral equation 

(40) MO = [ M£E[x]T ai )MP[^ Q1 )dT. 

Jo 

This shows that hoo is the fixed point of the particular smoothing transformation: 

T{hoo , /ioo) = Law (e[X] U ai V' + E[Y] U ai V") 


where U,V',V' are independent, U is uniformly distributed on (0,1), and V',V" are distributed 
with law hoo. Note that 


E[E[X]U ai +E[Y}U ai ] = E [ A 1 + E W = i 

1 + «i 

Positive fixed points of smoothing transformations are deeply studied, mainly in connection with 
the so-called Branching Random Walk, see e.g. [1, 37]. It is worth noticing that in general — with 
the exception of few cases, see [5, 6] — there is no analytical expression of the law the of the fixed 
point of a given smoothing transformation. Even in our case it does not seem easy to obtain an 
explicit expression of hoo■ An exception is the very special case when E[X] = E[Y} = 1, where it 
is immediate to see that the solution is given by the exponential distribution 


hoo(0 


1 


Nevertheless, it is always possible to recursively determine the exact expression of the integer 
moments of hoo- Indeed, using the fact that hoo = T(h. 0 0 , hoo) and the Newton binomial formula, 
if m-i := f v q hoo(dv) < +00 and toi = f vhoo(dv) = 1 , one gets 


(41) 


{E[X} + E[Y] - l)i + 1 - E[XY - E[Y] 


i -1 

E 

1=1 


ElXyElYy-^mjmi-j. 


As we shall see in Proposition 7 below, J v r hoo(dv) < +00 if and only if a r /r < a\. Hence, 
to* < +00 if and only if (E[X] + E[Y] — l)i + 1 — E[X] 4 — E[Y]* > 0 so that (41) is well-defined. 
Note that when at least one of the two moments E[X] and E[Y] is strictly bigger than 1, then there 
is r < +00 such that (E[X] + E[Y] — l)f + 1 — E[X] r — E[Y] r < 0 and hence / v r hoo{dv) = + 00 . 
This means that in this case h 0 0 has heavy tails. On the contrary if both E[X] and E[Y] are 
smaller than 1, then hoo Las moments of all orders. 

The fact that the steady state of (36) satisfies (40) suggests that, for a better understanding of 
the asymptotic behavior of equation (36), one can first to consider the large-time behavior of the 
equation 

(42) |M) + M) = /b(E[X]£)/i t (E[Y]£)., 

which is satisfied by the scaled solution h(£,t) = /t(e~ ait (/mo) of the equation 

(43) |/t( 0 + m = ft (e[A]£) ft (E[Y]£). 

Equation (43) is associated to the very simple linear collision rule 

(44) V' = YE[X] + VjE [Y], V- = V : E[X] + V 7 E[Y] 

and it is a special case of a general class of kinetic equations which have been deeply investigated, 
see for instance in [2, 3, 14, 40, 44], It is immediate to realize that the interactions (44) correspond 
to a mean version of the interactions (1), and can be formally obtained from (1) by substituting 
X and Y with their mean values. Of course, this is consistent only if the mean values of X and Y 
are integer numbers. For this equation, one can resort to the results of [2] to obtain 
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Proposition 7 ([2]). Let E[X + Y] > 1, and let (37) hold true for some r > 1. Then, there 
exists a unique solution to (40) such that +oq ~ ) vh oc (dv) = 1. Moreover, this solution satisfies 
fro+oo) v9 h°o(d v ) < +oo f or Q > 1 if and only if a q /q < a\. Finally, if Jvfo(dv) = mo < +oo 

and ft is a solution of (43), then f vft{dv) = moe ait and ft(e~ ait ^/m o) converges to /ioo(£) for 
every £ > 0 as t —> oo. 

Proof. The result follows from Proposition 2.1 and Theorem 2.2 in [2]. We have only to note that 
the spectral function fi(-) defined in Section 2.2 of [2] in this case coincides with r K»• a r /r — a.\. □ 

One can adapt some results in [2, 40] to obtain 

Proposition 8. Let E[X + P] > 1 and assume that (37) holds true for some r £ (1,2). Let f t be 
a solution of (43) with initial condition fo with finite r moment. Then 

(45) M r (f t ) < C^. 

Moreover, if h t (£) = /t(e -ait £/m 0 ), mo = / vfo(dv), and ^(^ is the same as in Proposition 7, 
then 

d* r {h t ,h 00 )<C 2 e^-™At, 

for a suitable constant C 2 = ^(/cu^oo)- 

Proof. The first part of the thesis is a minor modification of Theorem 3.2 in [40], observing that in 
this case the function 6(r) -defined in equation (32) of [40] - coincides with E[X] r +E[F] r —1 = a r . 
The last part of the thesis follows by a simple modification of Thm. 2.2 in [40]. See also Proposition 
3.12 in [4] □ 

In the next proposition we show that, under suitable assumptions, the scaled solutions of (2) 
and (43) merge when t -A + 00 . Combining this result with the known asymptotics for the scaled 
solution of (43) one obtains the asymptotic behavior of the scaled solution of (2) in Case 2. 

Proposition 9 (Scaled solution). Let E[X + Y] > 1 and assume that (37) holds true for 
some r £ (1,2). Let f^ be a solution of (21) with initial conditions fo with finite r-th moment 
and f^ be a solution of (43) with the same initial condition fo with mean mo- Set = 

ft l \e~ ait ^/mo) for i = 1,2. Then, 

4(/4 1} ,/4 2) ) < C(/ 0 )[te (Q '- r “ l)t + e-“ lrt ]. 

In particular, h{ 1 ' 1 converges weakly to the probability distribution hoo given in Proposition 7. 
Proof. We have 

r r |/t (1) (- mo)/P(- mo) - /t (2) (- E[x]c)/ t (2) (- e [Y]^j 
< r r |/ t (1) (- mo)/? } (- mo) - / t (2) (- mo)/? } (- mo) 

+r r |/t (2) (- mo )/* 2 } (- mo) - ft 2) (- Mo)/t (2) (w) | 

+ r r |/f } ( - kx(0)fi 2) (EMC) - / f (2) (E[X]c)/ t (2) (E[F]e) 

=: A\ + A2 + A3. 

Now, arguing as in the proof of Theorem 2, 

Moreover, using (68) 

a 2 < r r |/f } (- mo) ^ ft 2) ( e mc) | < c{Y)[i +M r (/ t (2) )] 

and analogously 

A 3 <C(X)[l + M r (f^)\. 
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At this stage recall that, by Proposition 8, 

M r (/ t (2) ) < 

Putting all estimates together one gets 


(46) 


m) ■■= r 


Now, consider that 

and 


fV (- mo) / f (1) (- mo) - ft 2) (- ms) fi 2) (- Erne 

< A r 4(/ t (1) ,/ f (2) ) + C 3 [l + e«i. 

|(e7f } (0) = e‘/ t (1) ( - k x (0)fT\- k Y { 0), 
|(e7t (2) (0) = e7 t (2) (E[A]c)/ f (2) (E[F]c). 


Last, consider that /q 1 ' = /q 2 \ and, using (46) 


Jr 


/?>«) - /t 2) (0 < f e s H s (£)ds < f e a [A rl r r (fP,fW) + C 3 (l + e“' 8 )]cfo. 

Jo Jo 


Hence, 


e*4(/ t (1) ,/ t (2) )<A r / + / [ e * + e< 1+ “'>»)]ds. 

Jo Jo 

By Lemma 18, 

4(/t (1) ,/ f (2) ) < c r [M r (/ t (1} ) + M r (/ t (2) )] 

and, by Propositions 5 and 8, t i—>• M r (f 7) + M r (fj 2> ) is a locally bounded function. Hence, 
Gronwall Lemma gives 


<(/ t (1) ,/ f (2) ) < /V“ Ar)s + l)]ds. 

JO 


If now a r ^ 0, it holds 


d* r {h ( t 1) ,h ( t 2) ) < C 3 e~ airt [te art + 1 ——] < C 4 [le(-“ ir+ “’-) t + e- airt ]. 


In the remaining case a r = 0 it holds d*(h[ 1 \h[ 2 ' > ) < C 3 e airt . 


□ 


5. The quasi-invariant limit 

In this section we give precise statements and proofs of the limit procedure considered in Section 
2.2. 

Proposition 10. Let ft,e(z) be the solution of (4) for X and Y defined in (8) where E[A 2 +Y 2 ] < 
+oo and /o has finite second moment and mean equal to m 3 . Then, g(t,z) = lim e _>.o ft/e,e( z ) 
satisfies 

d t g(t, z) = bx(Px(z) - z)d z g(t, z) + b 2 m. 0 e ait (py(z) - 1 )g(t, z ), 

(47) 

<?(M) = f 0 (z). 

for ( t , z) £ (0, +oo) x [0,1], where a.\ = &i(E[X] — 1) + & 2 E[F]. Moreover, the initial value problem 
for equation (47) has a unique solution, and 

Mi(g(t)) = m 0 e ait . 

Proof. Existence and uniqueness of a unique bounded solution in Q = (0,1) x (0 ,T), for every 
T > 0 to equation (47) is known (cf. [25] for a detailed study of this type of linear equations). 
Simple computations show that 

a 2 := E[X] 2 + E[P] 2 - 1 = 2e(E[X] - l)h + e 2 ((E[A - 1]) 2 & 2 + (E[+]) 2 ^), 

while 

P := Var(X) + Var{Y) = e(E[(l - 1) 2 ]&! + E[F 2 ]6 2 ) - e 2 ((E[l - 1]) V + (E[F])V), 
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and 

7 := E[X]E[F] = e& 2 E[f](l + e 6 i(E[X] - 1)). 

Recall that a\ = E[X] + E[T] — 1 = ear. Hence, owing to (32), one finds that 

(48) Var(f t ,e) < G(et), 

for a suitable function G bounded on any fixed interval of time [0,T]. In addition we have 

\d 2 JtA Z )\ = l)z V ~' 2 ftA V )\ < ^2 v2 ft,e( v ) = Var (ft,e) + Mi (/ t , e ) 2 • 

v>2 v>2 

Recalling that that Mi(/^ e ) = moe aiet and using (48), for every T < +oo one gets the bound 

sup \d 2 Jt, e (z)\ < c, 

t:te<T 

where C = C(/ 0 , X,Y,T) is independent of e < eo and z. Hence d z f t , e {z) is locally uniformly (in 
t and e) Lipshitz in z, which means that 

(49) d z f t<e (z)\ z = Zo - d z f t:e (z)\ z = Zl < L(T)\z 0 - zx\. 

for every zq, Z\ in [0, 1 ] and every t > 0 such that te < T and e < eo- 
To prove convergence, write 

ft(Px{z)) = f t {z) + ebi(j)x{z) - z)d z ft(z) + eb 1 R lje (t,z), 

and 

ft(pr(z)) = f t ( 1 ) + eb 2 (py(z) - l)d z f t (l) + e& 2 R 2 j£ (f, z). 

In the previous equations, 

Rl,e(t,z) = ( Px(z) - z)(d z ft(Zl,e,t) ~ 9 z f t {z)), 
for some 2 i j£jt £ [ 0 , 1 ] such that |z — < ebi\p-^(z) — z | and 

^ 2 , e (t, z) = (py(z) - 1 )(d z ft(z 2 ,e,t) ~ d z f t { 1)), 
for some Z 2 , e ,t such that 0 < 1 + b- 2 e(py(z) — 1) < z 2 , £jt < 1. Consequently 

ft(Px(z))ft(P Y {z)) - f t (z) = eb 2 f t {z)(py(z) - 1 )m{t) + eb 1 (p ji (z) - z)d z f t (z) + ei? 3 , £ (f, z). 
Here, m(t) = d z f t { 1) and 

R 3 , e (t,z) = ebib 2 (py(z) - z)(py(z) - l)d z f t (z)m(t) + ebib 2 {p^{z) - z)d z f t (z)R 2 , e (t,z) 

+ b 2 R 2 ,e{t , z)ft(z) + Ri,e{t, z)bi[l + eb 2 (py(z) - 1 )m(t) + eb 2 R 2 ^(t,z)). 

Let us set g t ,e(z) = ft/e( z )- Then 

(50) d t g t , e (z) = g t ,e(z)b 2 (py(z) - l)Mi(g ti£ ) + bi(px(z) - z)d z g t , e {z) + R 3 )£ (f / e , z )- 
Using (49), for a suitable constant C'i(T) we obtain 

sup Y |-Ri, e (V e ^)| < eCx(T), 

(z,t)eQ i=12 

so that 

(51) sup \R 3 , e (t/e,z)\ < eC 2 (T). 

(z,t)£Q 

Recalling now that for every z £ [0,1] and every t £ [0, T\ 

d z gt,e ( Z ) < M l(gt,e ) < C 3 (T), 

using this with (51) in equation (50) we obtain \d t g t ^{z)\ < C 3 (T). Summarizing, for a suitable 
bounded positive constant C${T) we have 

sup \d z g t , € {z)\ + sup \d t g t ,e{z)\ < C 5 (T). 

(z,t)e[0,l]x[0,T] (z,i)e[0,l]x[0,T] 
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Combining Ascoli-Arzela theorem with the Banach-Alaoglu-Bourbaki theorem one concludes that 
there is a subsequence g., €n (') converging strongly in C°([0,1] x [0, T]) (for every T) and weakly-* 
in W 1,00 (Q) to a function g. 

Since R^ e —» 0 uniformly by (51), one can pass to the limit in (50) to obtain that g is a solution 
of (47). Next, uniqueness of the solution of (47) gives that the whole family g. e converges to g 
strongly in C°([0,1] x [0,T]), for every T > 0. Finally, since the mean and variance of g t<e are 
uniformly bounded with respect to e, the sequence {gt,e}e>o is tight, and the limit function g(-,t) 
is also a p.g.f. with finite variance. □ 

Define 4>x{z) = 12k>o z k P{X > k} and <j>y(z) = J2k> 0 z k P{Y > k}. Setting to = E[X], if 
&i(E[X] — 1) + 6 2 E[T] = 0, one can write (1 — 771 ) 61/62 = E[F] so that P*(z) := m~ l (j>x{z) and 
Q*(z) := 62 / 61(1 — m)~ 1 <j)Y(z) are p.g.f.. It is worth noticing that, if b 1 = b 2 = 1, P* and Q* are 
the p.g.f. of the so called renewal distribution of px and py. Recall also that given a probability 
density g on N with finite mean M\(g) = kg{k), its size-biased version is defined by 


9*{k) 


kg(k) 

Wa)' 


Proposition 11. Under the same assumptions of Proposition 10, i/ 6 i(E[A] — 1) + & 2 E[y] = 0, 
then the unique C 1 -function goo{z) with 9 oo(l) = 1 which is a stationary solution of (47), i.e. the 
unique solution of (16), is given by 


(52) 



(1 -Py(s)) , 
(Px(s)-s) 


}■ 


Moreover, g^ is the p.g.f. of an infinite divisible distribution g oo over the non-negative integers. 
Finally, goo is the density of a random variable V/o satisfying 

(53) C^ + l + C + f^fc 

k =1 


where is the size biased version of V/o, f has p.g.f Q*, fk has p.g.f P* for every k > 1 , 
7 has geometric distribution of parameter m = E[A] and V/o C;7i Cu ^ 2 , • • • ore stochastically 
independent. 


Proof. A stationary distribution of (47) must satisfy (16). That is 

/n m 0 6 2 (l -Py{z)) ^ 

dz9oo{z) = , , - r~9co\ z ), 

bi{Px(z)-z) 

that gives (52). Now note that, using (24) and the analogous relation for Y, one gets 


R{z) 


6 2 (1 — Py(z)) _ b 2 (t>Y{z) 

bi{Px{z)~z) 6 i(l - (fx(z)) 


(1 - m)Q*(z) 
1 - mP*(z) 


G m {P* (z))Q* (z) 


where 


is the p.g.f. of a geometric random 
the random variable 


G m (s) = 


1 — TO 


1 — ms 

variable of parameter to. In other words R(z) is the p.g.f. of 




k =1 


We proved that 


9oo(z) = exp | - m 0 J R(s)ds|, 


where R is a p.g.f. and hence m$R is an absolutely monotone function. This show that goo{z) is 
a p.g.f. of an infinite divisible distribution goo over the non-negative integers. Indeed a function h 
with h( 0 ) > 0 is the p.g.f. of an infinitely divisible distribution over the on-negative integers if and 
only if h(z) = exp{J/ H(s)ds} for an absolutely monotone function H (cf. Thm 4.2 [45]). Finally, 
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note that m 0 1 d z g 00 {z) = z 1 g% 0 {z) with g^z) = Efc>o zk 9to( k ) and 9lo( k ) = kg°o(k)/m 0 , that 
is the size-biased version of g^. At this stage one can re-write (16) as 

zR{z)g 00 {z) =g* 0 {z). 

Therefore, in terms of random variables, (16) takes the form 


(54) 


C = Ko + i + C + 


k =1 


where is the size biased version of V^. 

Adapting Lemma 2.2 in [18] one gets a further result. 


□ 


Lemma 12. Let ( f n ) n , (<?n)n>/i<? be the Laplace transform of probability measures on the positive 
real axis, such that f n —> /, g n —>• g, Mi(fn) = M\{g n ) = mo, and for some e > 0 


(55) 

Then 

(56) 


sup | max | f v r+e df n , f u I ’ +€ dg ri }} = M e < +oo. 

n 1 1 Jr+ J r + ' > > 


d*{f,g) < liminf d*(f n ,g n ). 


Proof. Recall that pointwise convergence of Laplace transforms yields uniform convergence on 
every compact set, hence, for every 0 < 5 < R < +oo, 


(57) 


c„ = C n (S, e) := sup |/„(£) - f{f)\ + sup \g n {£) - g(£)\ 0 . 

£e[<5,K] £e[<5,-R] 


Moreover, since pointwise convergence of Laplace transform yields weak convergence of distribution 
(cf. Theorem 25.11 and Corollary 25.12 in [12]), using condition (55) one obtains Mi(f) = Mi(g) = 
m.Q and 


max 


{ j v r+e f(dv),J v r +‘g(dv)}<M e . 


Using again (55), (64) yields 

i mo - m\r r < i mo -1 - rn. 0 f\r r + \m - 1 - m 0 m~ r < 2c, +e m 

Hence 

sup |/„( 0 -/( 0 ir r <2c r+e 5 e M e . 

«e(o,<5) 

Now, consider that 

d*(f n ,f)< max { sup | M0~f(0\C r , sup \f n (0 - /(£)|£ _r , sup \f n (0 - /(OI£ _r j 

L «e(o,5) «e[<5,R) £>R J 

< max 12 c r+e 6 € M e , C n , 

Since the same estimate holds for d*(g n ,g), 

lim [d*(f n ,f) + d*(g n ,g)] =0. 

n 

To conclude it suffices to observe that 

d*(f,g) < d*(f n ,g n ) + <(/„,/) +d*{g n ,g). 


□ 


Theorem 13. Let X and Y be defined as in ( 8 ), with bi = & 2 > a nd X and Y satisfying (7)-(31). 
Let f^l, i = 1,2, be two solutions of (2) with initial conditions f^ with mean equal to mo and 

finite second moment. Let g ^ and g^ be the corresponding quasi invariant collision limits of 
Proposition 10. Then, for every 1 < r < 2, 

<( 5t (1 U (2) ) < d;(g£\gW)e-^ [ x ]r -H' 
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Proof. We first observe that E[X] + E[Y] = 1 implies 1 — E[X] > 0. Explicit computations and 
Taylor formula give 

a r = E[X] r + E[Y] r - 1 = [1 + e(E[X] - l)] r + e r E[X] r - 1 


(58) 


= re(E[X] - l) r_1 + 


-1 i r ( r 1) ^irpr id i /1 i \r —2 , 


2 (E[X] - 1) 2 (1 + Qff- 1 + e r E[X] r 


with \Q e \ < e(l — E[X]). Therefore 

(59) a r = re(l — E[X]) r + o(e). 

Combining (35), (5) and (59) it follows 

M 2 (f^ et )<C[l + e- a 2 ^]<C', 

for every t> 0 and e > 0. Since g[ l \z) = lim e ^ 0 ft/\ € (z), combining Theorem 2, (56) and (59) 
«C(St (1) >flt 2) ) < liminf < d*{g£ ] ,^ 2) )e“ r(1 " E[x])r 1( . 


□ 
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Appendix A 

A.l. Recurrence relations. We state without proof two previously used results. 

Lemma 14. Let 0 < ao < +oo, and, for every n > 1, let a n > 0 satisfy, for some 0 < A < +oo, 
the inequality 

A ”” 1 

(60) a n < — ^ aj. 

71 j= o 


Then, for every n > 1 
(61) 


a n < a 0 - 


T(A + n) 


T(A)r(n+l)' 

Moreover the equality sign in (60) implies the equality sign in (61). 


Lemma 15. Let 0 < ao < +oo, and, for every n > 1, let a n > 0 satisfy the inequality 


A 


n—1 


< , n H 


J=0 

for some 0 < A, B, /3 < +oo. TTien there is K = K(A, B, (3, ao) > 1 such that for every n > 0 

a n < K n+1 . 

We conclude this short section by recalling that, for any t > 0 and A > 0 [30] 


(62) 


= (A-l)t _ 


= ^e- t (l-e- t r 


n>0 


r(A + n) 
r(A)r(n + i)' 


A.2. Inequalities for Laplace transforms and cumulants. Given a random variable Z taking 
values in R + let us set, for £ > 0 

Lz(0= K [ e ~ iZ ] and k z (Q = log(L z (0)- 


Since for all x > 0 

\e~ x - 1 | < M, 

and for every r € (1, 2) one can find a constant c r such that 

(63) \e~ x — 1 + x\ < c r \x\ r , 
whenever E[Z r ] < +oo and £ > 0 it follows 

(64) IMO-l-mi <c r enz r ], 
and, if E [Z] < +oo, 

(65) \L Z (0 - 1| < 4E[Z] 

The following lemmas contain useful bounds for cumulants. 


Lemma 16. Let E[Z r ] < +oo for some r £ (1,2). Then, if 0 < £ < (2E[Z[) 1 , 

( 66 ) \k z (o+m[z}\ < c r f r E[z r }+e(nz}) 2 . 

Moreover 

(67) S u P ^#^=E[^] r . 

C>o Isl 

Proof. For every complex number 3 write 

log(l + «) =z{l + R{z)). 

Then |1?(^)| < \z\ for every 2 such that \z\ <1/2 (cf. Proposition 8.46 in [15]). Thanks to (65), 
0 < £ < (2E[Z']) _1 implies 


\Lz(o-M<m[z\<l- 
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By virtue of (64) and (65) 

M£) + m\ = Ilog(l + L z (0 - 1 ) + 0E[Z]| = I (L z (0 - 1)(1 + R(L z (t) - 1)) + &[Z}\ 

< i Lz(o - 1+ sew + [l z (o -1] 2 < x r cnz r }+ emf- 


To obtain (67), note that £ i—>• L z (£f) = E[e ^ z ] is a non-increasing function and E[e ^ z ] < 1. 
Hence £ i-A | log(E[e _?z ])| = \kx{£,)\ is non-increasing. So 


sup 

t>o 


jMor 


lim 


jMor 

ki r 


Thanks to the previous estimate of the Lemma 

Mor . mz]+pjoi 

|^|r |£|r 

with |p(£)| r < C|£| 2r . This gives (67). □ 


Lemma 17. Let Z be a random variables taking values in N and X a non-negative real random 
variable. Assume that E[X r ] +E[i? r ] < +oo for some r £ (1, 2). Then there is a positive constant 
C = C(X) such that, for all £ > 0 

( 68 ) | L z (-kx(0) ~~ L z mmc r < C{X)[ 1 + E[Z r ]. 

Proof. Let us define 

C(X) = max{c r (E[X r ] + 2E[X] r ) + 2 r ” 2 E[X] r , 2 r+ 1 E[X] r }. 

If £ > 1/(2E[X]), it holds 

\L z (-k x (0) - L z (fE[X})\r r < 2 r+ 1 E[Xf < C(X). 

If now £ < 1/(2E[X]) one obtains 

\L z (-k x (z))-L Z mmc r < i L Z (-kx(Q) -1 - k x (mz]\r r 

+ |1 - ME[Z] - Lz(E[X]Q\C r + E[Z]mX] + k x m~ r 

< CrE[Z r ]{C r \k x (Z)\ r + E[X] ? '} + E[Z r }{c r E[X r ] + £ 2 ~ r E[X] 2 } 

< E[Z r ][c r {E[X r ] + 2E[X] r ) + 2 r " 2 E[X] r ] 

Note that the last two inequalities from below have been obtained by using (64) and (66) and 
(67), respectively. □ 

Lemma 18. Let C := sup ?>0 |1 — e -? ||£| -1 < + 00 , and B = sup 0 < j< lo g( 2 ) |£||1 — e _? | _1 . Then, 
for every r > 1 

d* r (f,g)<C r d r (f,g), 

and 

d r (f,g) < max{2 r+1 , B r d*(f, g)}. 

Moreover, if M\{f) = M\ (g), and M r (f) and M r (g) are finite for some r £ (1,2], then 

d*(f,g) < Cr[M r (f) + M r (g)], 

for a suitable constant c r . 

Proof. Noticing that /(e - ^) = /(£), one can write 

l/(£) ~g(£)l = l/(e~')~g(e~^)| 11 — e" g l r 

ki r ii-e-«i r kr 
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Hence d*(f,g) < C r d r (f,g). Moreover, 


\f(s)~g(s)\ 


d r (f,g)=m ax< sup — 

I s6(0,l/2) I 1- 


, sup 

«e(i/2,i) |1 s l r 


< max l 2 r+1 


|/(e *)-g(e «)| |£| r 


€6(0,log(2)) 


Il-e-€| 


< max { 2 r+1 , B r d* (/, g )} . 


Thanks to (63) 


\f(0~g(0\ = | f(e~ iv - 1 + vQf(dv) - f(e~t v - 1 + Qg(cfa)| 
|^|r |e 

< c r [M r (f) + M r (g)\. 


A.3. Moments of random sums. Let V\ and Vi two integer valued random variables. Let 
(Xk)k> l and (Y k )k> i two sequences of iid non-negative random variables and assume that V\. Vi, 
{X k )k>i,(Y k ) k >i are independent. 

Lemma 19. Under the previous assumptions let r £ (1,2] and assume that E[y r + Y r ] < +oo 
and E \V{ + Vf] < + 00 . Then, if r = 2 
, Vi u 2 „ 


E[( g Xi + g g ] = E[H 1 ]Far(X 1 ) + E [V^Var^) 


i= 1 2 =1 


+ 2E[H 1 ]E[y 2 ]E[X 1 ]E[y 2 ] + E[H 1 2 ]E[A ' 1 ] 2 + E[y 2 2 ]E[y !] 2 


while if r £ ( 1 , 2 ) 


Vi v 2 


E[ (g x *+g y,) r ] <e[h 1 ]e[a[] + e[h 2 ]ek'-; 


2=1 2=1 


( 2 E[y 1 ]E[H 2 ]E[X 1 ]E[y 1 ]) 3 + E[y i r ]E[A 1 ] r + E[y 2 r ]E[y 2 ] r . 


Proof. We only proof (70). Let rq and n 2 two integer numbers. Since r/2 < 1, for any given set 
of positive numbers z \,..., z n , (X^=i 2*g 2 — 1 Z J 2 ■ This implies 


711 Tl<2 OH. ' fc l 2 — , _ j_ 

®i(E**+E«) VC 1 ) +(Et) +2E^Et.)’i 

2=1 2=1 2=1 2=1 2=1 2=1 

ni n .2 ni 222 ni n 2 r 

= E[(gx 2 + gy i 2 + 2 gx i gy i + g g yg) 5 ] 

i—1 i=l i—1 i=l i^j,i,j=l i^j t ij =1 

ni 222 ni 722 r ni r 712 

<E[gx[ + gyr + 2 s(gx i gy i ) 3 + ( g M,g + ( g 


72 1 n 2 


2=1 2=1 

ni 222 


221 222 


2=1 2=1 2=1 2=1 iy£j,i,j = 1 




221 222 


<E[^X[ + ^y; + 2 i 


2=1 2=1 


2=1 2=1 


E 

* Am, 1=1 


E y < y j 

*A‘,*iJ= 1 


< mE[A[] + n 2 E[y i r ] + 2 $E[^ X, ^ E;] 5 + E[ ^ X ?; X,]I + E[ E 

*=i i=1 *A'i*iJ= 1 *A',*,j'= 1 


< niE[X[] + n 2 E[y 1 r ] + 2? (nm 2 E[X 1 ]E[Yi]) 3 + n , 1 'E[X 1 ] r + ?r^E[y 1 ] r . 

Note that we used Jensen’s inequality to get the last line. To conclude let us apply the previous 
inequality conditionally on V) and y 2 to get 

Vi v 2 r 

E[( E ^ + g e) ' ] < E[Vi]E[Xf] + E[Vi]E[Y{} + E^JE]^ 5 ] ( 2 E[X 1 ]E[y 1 ]) 3 + 


2=1 2=1 


E[y 1 r ]E[x 1 ] r + E[^]E[y 2 ] r 
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and then apply the Jensen’s inequality one more time to get Efl^ 2 ] < E[V)] 2 . □ 

A.4. Gronwall’s inequality. A well-known version of the Gronwall lemma is the following. 


Lemma 20. Suppose that 1 1 —> u{t) and t i-a h(t) are integrable functions on the interval [0,T]. If 
for some constant C > 0 


u(t) < h(t) + C f u(s)ds 0 < t < T, 
Jo 


then, for every t £ [0, T] 


u(t) < h(t) + C 



A simple consequence of the previous Lemma is the following. 


Lemma 21. Let t <—> u(t) be a positive measurable function on the interval [0,+ 00 ) such that: 

(i) I = {t : u(t) < + 00 } is an open interval (0, so), 

(ii) u is continuous on I 

(iii) if so < +00 then 

lim u(s ) = + 00 . 
t->s 0 

If there is a non-negative continuous function t 1 —> h(t) on [0, + 00 ) and a constant C > 0 such 
that 

u(t) < h(t) +C u(s)ds , 

Jo 

for every t > 0, then 

u(t) < h(t) + C f e C{ - t ~ s ' , h{s)ds, 

Jo 

for every t > 0 and, in particular, sq = 00 . Finally, if h(t) is of bounded variation then 

u(t) < e Ct h(0) + f e c( ' t ~ s ^dh(s). 

Jo 

Proof. It is enough to prove that So = + 00 . Indeed in this case the thesis follows from the previous 
Lemma. Assume that so < + 00 . Let H(t) = h(t) + C fg e c ^~ s ^h(s)ds. By (ii) and the fact that 
H (t) is bounded on every finite interval, one can choose e in such a way that u(so — e) > H (so — e). 
Now u is continuous in (0, so) and hence it is bounded on (0, Sq — e). Hence, by the previous 
Lemma, u(so — e) < H{so — e) which gives a contradiction. The last part of the proof follows 
resorting to integration by parts. □ 
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